Perfect magnetohydrodynamics as a field theory 
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We propose the generally covariant action for the theory of a self-coupled complex scalar field and 
electromagnetism which by virtue of constraints is equivalent, in the regime of long wavelengths, to 
perfect magnetohydrodynamics (MHD). We recover from it the Euler equation with Lorentz force, 
and the thermodynamic relations for a prefect fluid. The equation of state of the latter is related to 
the scalar field's self potential. We introduce 1 + 3 notation to elucidate the relation between MHD 
and field variables. In our approach the requirement that the scalar field be single valued leads to 
the quantization of a certain circulation in steps of h; this feature leads, in the classical limit, to 
the conservation of that circulation. The circulation is identical to that in Oron's generalization 
of Kelvin's circulation theorem to perfect MHD; we here characterize the new conserved helicity 
associated with it. We also demonstrate the existence for MHD of two Bernoulli-like theorems for 
each spacetime symmetry of the flow and geometry; one of these is pertinent to suitably defined 
potential flow. We exhibit the conserved quantities explicitly in the case that two symmetries are 
simultaneously present, and give examples. Also in this case we exhibit a new conserved MHD 
circulation distinct from Oron's, and provide an example. 

PACS numbers: 95.30.Qd, 52.65.Kj, 47.75.+f, 47.65. +a, 03.50.-z 



I. INTRODUCTION 

It is well known that perfect fluid dynamics admits a 
description in terms of dynamics of a self-interacting real 
scalar field. Analogously, a superconductor, insofar as 
we are interested in the flow of current through it, can 
be viewed as a complex (charged) scalar field with self- 
interaction coupled to the electromagnetic field yj. Can 
the flow of a plasma interacting with the electromagnetic 
field be treated analogously ? In one approach such a 
plasma is represented by two charged Schrodinger-like 
scalar fields representing the flows of ions and electrons, 
respectively 0- 

Here we ask, is perfect magnetohydrodynamic (hence- 
forth MHD) flow, namely flow of plasma with negligible 
charge separation, and with a magnetic field frozen into 
the flow, amenable to relativistic description in terms of 
the dynamics of a single, possibly complex, scalar field 
interacting with the electromagnetic field ? Such a de- 
scription would not only be of methodological interest, 
but might supply new insights into MHD flow as well as 
hints about solution of its intricate equations in a wider 
range of problems than possible hitherto. 

In this paper we approach the problem by writing a 
joint action for the electromagnetic field and a complex 
scalar field, the latter self-coupled as well as coupled to 
a non gauge vector potential. Part of the action, which 
is generally covariant, involves constraints, one of which 
enforces the condition for the magnetic field to be frozen 
into the flow. We show that in the long wavelength limit 
the dynamics following from this action is that of MHD. 
Among the immediate consequences of the field point of 
view are concise derivations of the circulation conserva- 
tion theorem in MHD, and of Bernoulli- like theorems, 
one for each spacetime symmetry. 

Our approach offers a subtle resolution to a quandary 



originating in scalar field formulations of hydrodynamics. 
Let us understand the problem. The energy-momentum 
tensor for a relativistic real scalar field whose gradient is 
timclike can be cast in the form of the energy-momentum 
tensor for a fluid whose density and pressure are related 
to the self-interaction potential, and whose 4-velocity is 
the gradient of a function of the scalar field. Thus only 
potential flow is so described. It would seem that a scalar 
field representation of hydrodynamics is incapable of cov- 
ering all types of flow, even if we grant the perfect flow 
condition. 

It is interesting to contrast this drawback with a simi- 
lar one evidenced by the purely fluid Lagrangian formu- 
lation of isentropic hydrodynamic flow [310, 0J IE The 
velocity field is there the sum of two terms: the gradient 
of the Lagrange multiplier responsible for enforcing mass 
conservation, and the entropy per unit mass times the 
gradient of the Lagrange multiplier responsible for en- 
forcing entropy conservation (which is a prerequisite for 
isentropic flow) . The velocity field is thus again potential. 

The above problem is usually fixed by invoking Lin's 
trick 0: introduce a new Lagrange constraint enforcing 
the conservation of some additional local quantity along 
the flow. In this case the velocity field is generic and 
not restricted to potential flow However, the 

nature of the extra conservation law is mysterious. It 
is usually claimed that the conserved quantity is one of 
the Lagrangian coordinates of fluid elements. But then 
the question arises, why only one of the three such coor- 
dinates enters. And even if there were a principle that 
chose one of the three for its special role, one may ask, 
why cannot the problem be formulated entirely in Eule- 
rian coordinates ? 

In the approach proposed here the fluid 4-velocity is a 
scalar quantity times the difference of the scalar field's 
gradient and the 4-potential (not a gauge one). The 4- 



2 



velocity is not collinear with any gradient because the 
phase, being a modular variable, does not necessarily 
yield a true gradient (one whose curl vanishes) upon dif- 
ferentiation. This puts both vortical and potential flows 
within the province of the theory. 

In Sec.[n]we define the scalar and electrodynamic vari- 
ables of the theory, and propose the action including the 
constraints which establish the connection with MHD 
flow. In Sec. IIIII we derive from the least action principle 
the theory's equations, including the source of Maxwell's 
equations and the MHD Euler equation of motion. We 
also discuss the gauge freedom of the theory, showing that 
it is limited to a U(l) x U(l) group. In Sec. II VI we obtain 
the equation of state of the fluid represented by the scalar 
field and its fluid-like energy-momentum tensor. Sec. 
provides a translation of the theory into 1 + 3 language 
in order to clarify the theory's content. For example, we 
show how to represent the kinematical variables of the 
flow in terms of the scalar's phase. Sec. lVIl is devoted to 
Oron's generalization of Kelvin's circulation theorem to 
perfect MHD [HG3 in the scalar field language; it also 
characterizes the new conserved helicity associated with 
it. Finally, in Sec. IVIll we demonstrate the existence in 
MHD of two Bernoulli-like theorems for each spacetime 
symmetry of the flow, one of these applicable to suit- 
ably defined potential How. We exhibit the conserved 
quantities explicitly in the case that two symmetries are 
simultaneously present, and provide examples. Sec. IVIII 
also exhibits a new conserved circulation distinct from 
Oron's and Kelvin's, and provides an example. 

We work in units with c = 1. Our signature is 
{ — ,+,+,+}. Symmetrization and anti-symmetrization 
of tensors is denoted by Tr a p> = ^(T a + Tp a ) and 
Tj a/ 3] = \{T a p — T(j a ), respectively. In like manner 
W[ a/ 3-j] = \{W afs . n + W ia .f3 + Wp T , a - (a (3)}. 

II. THE THEORY'S CONTENT 
A. Fluid 

The four velocity of the MHD fluid is written u a while 
its baryon proper density will be denoted n. Other ther- 
modynamic quantities are the energy proper density p, 
the (isotropic) pressure p and the relativistic enthalpy 
per baryon fj,, 

B. Electromagnetic field 

Relatively to u a , the Faraday tensor F a p and its dual 
*F a f3 can always be decomposed into electric and mag- 
netic parts, 

F a p = u a Ef3-Uf 3 E a +e a p 1 su l B 5 , (2) 
*F a p = -u a Bi3 + U{3B a +E a p 1 su l E s ', (3) 



where the electric field is defined as E a = F a p u" and 
the magnetic field is B a = ^E a 0j5 u^F^ 5 = *F 7Q u 7 , 
respectively. The Levi-Civita tensor is normalized such 
that e i23 = \F r 9- 

In perfect MHD the conductivity of the fluid is as- 
sumed so high that it completely suppresses the electric 
field: E a — 0. Thus in MHD the Faraday tensor is char- 
acterized by the condition 

F a0 u 13 = 0. (4) 

For performing variations the Faraday tensor must be 
expressed in terms of potentials. Our notation is F a p = 

A/3-a — A a .p — A/3 t a — A a fi. 

C. Scalar field 

We shall represent the fluid by a complex scalar field 
ip whose basic action is taken to be 

where rj a is an auxiliary vector field. We have added 
a self-interaction potential V(ipip*) in accordance with 
experience in representing fluids by scalar fields. Why is 
this an appropriate action ? Were ip here an electrically 
charged field (charge e), we would write its coupling to 
the gauge potential A a as above, but with r\ a i— > eA a . 
An MHD fluid is a mixture of two oppositely charged 
fluids, e.g. electrons and ions, so each component would 
deserve an action like the above, but with opposite signs 
of e. Since MHD assumes that the net charge in each 
small volume vanishes, we must contemplate a (perhaps) 
imperfect cancellation of the gauge interaction. Thus if 
we represent the MHD fluid by a single complex scalar 
field ip, this last must no longer be minimally coupled to 
A a . However, we know physically that ip must have a 
residual electromagnetic interaction. The most economic 
way to write its effective action is as in Eq. |JS), where 
f] a is expected to bear some complicated relation to the 
electromagnetic field [to be deduced in Eq. I|18|l below]. 
To put S^p in final form we write tp — a e lip (a and ip 
real) . In addition to the elements already mentioned, we 
shall have need of an additional vector field b a to play 
the role of a Lagrange multiplier charged with the task 
of enforcing the MHD condition |@J|. 

D. Action 

Omitting gravity's dynamics, the action of the theory 
is S — Sm + S$ + S c + S nor + where 

Su = ~ / F^Fcpi-gf/Wx, (6) 

Sv- = -- f[v, a cj ta + a 2 { V a -i 1 a ){^, a ~i la ) (7) 
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S c 



+V(a 2 )](-g) 1 / 2 d 4 x, 

^J(F aP + Kg aP )u a b^(-g)^ 2 d 4 x, (8) 

^J^n(u a u a + l)(-g)^ 2 d 4 x, (9) 

C / nu"i|„ (-g) 1/2 d 4 x. (10) 



5m is the usual Maxwell action; it is to be viewed as 
functional of the electromagnetic 4-potential A a . In ac- 
cordance with our earlier remarks, the scalar field action 
S^, now rewritten in terms of a and ip, represents the 
MHD fluid. We do not include a kinetic term for r\ a 
principally because, as we shall see, the apparent gauge 
invariance of the theory is a mirage. 

The constraint part of the action S c takes care of the 
MHD condition Indeed, interpreting b a as local La- 
grange multiplier, we see this condition emerges from 
variation of the first term in S c with respect to b a . But 
this is not correct in itself. By antisymmetry of F a p, 
F a p vP u a — 0, i.e. the vector F a p it' 3 has only three inde- 
pendent components. Hence the MHD condition should 
be derivable with help of a triplet of Lagrange multipliers; 
a 4-vector of them is too much. Accordingly we subject 
b a to the constraint g a p b a — by way of inclusion in 
S c of the term with local Lagrange multiplier re. Thus 
only three of 6"'s components are independent. 

The 4-velocity of the fluid u a is to be determined by 
the theory in terms of scalar field variables. We have 
no guarantee that it will be properly normalized, so we 
impose such normalization as a Lagrange constraint by 
the part of the action S nor . The £ is a local Lagrange 
multiplier and the factor n is included so that 4-velocity 
normalization is not required in regions not occupied by 
the fluid. 

Addition of to the action is found to be necessary 
for the correct relation between u a and the field variables 
[see Eq. <|22[l ] to arise; we have found no illuminating in- 
tuition for it. In Sq the C, is a constant with dimensions of 
action; it will emerge in Sec. lVI Al that C can be identified 
with the quantum of action H. 



III. EQUATIONS 

Variations will be carried out with respect to the inde- 
pendent fields: A a ,r] a , ip, er, 6 Q ,u a , g a p and n. We stress 
that for this purpose u a ,b a , A a and r\ a are to be varied 
independently of the other variables, but e.g. variation 
of u a entails a variation of u a and of g a p- 



A. Constraints and conditions 

Variation of S with respect to £ gives n (u a u a + 1) = 0; 
thus 



Accordingly in regions occupied by the fluid (n / 0), u a 
is normalized like any other 4-velocity. 

Variation of S with respect to re as well as b@ gives 



u a b a = 0; 
F a p u a + ku/3 = 0. 



(12) 
(13) 



Contracting Eq. (|13fl with u° gives nu^ up = 0; com- 
paring with Eq. i|llfl we conclude that wherever n ^ 0, 
k = 0. Thus Eq. I|13|) is equivalent to the MHD condition 
0. 

Variation of S with respect to u a gives 

^-F^b 13 -^nu a -C,nj la = Q. (14) 

47T 

Contracting this with u a and using Eqs. (@J and (f 1 1|) 
leaves us with 



ntt-Cv* u a ) = 0. 



(15) 



On the other hand, we may vary S with respect to n to 
obtain 

n[i((« a M a + l)+(i) a «1 = 0. (16) 

Comparison with Eq. i|15|) shows that wherever n/Owe 
must have £ = as a consequence of Eq. I|ll|) . Combin- 
ing with our previous result on re, we have for the local 
Lagrange multipliers 



k = £ = where n ^ 0. 
Now solving Eq. i|14fl for r\ a yields 

F a pbP 



47r£ n 



(17) 



(18) 



from which it is obvious that (we shall prove presently 
that ( ^ 0) 



i P =Vpb 13 = where n ^ 0. (19) 



u a u a = —1 where n^O. 



(11) 



One further constraint is obeyed by r^. By the MHD 
condition Q F a p as given by Eq. (0 can be cast as 

F a p = s a p lS u~< B s . (20) 

Using this form in Eq. i|18fl shows that 

rjp B p = where n ^ 0. (21) 

We see that the "residual gauge field" r\ a is heavily con- 
strained [Eqs. (HHJ and |JST}]. 

B. Equations for scalar field 

We shall now vary S with respect to r] a obtaining 

^^r 1 . 2 ^"-^); (22) 
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this equation links the hydrodynamical variables n, u a 
with the scalar field variables <p, a. If we now vary with 
respect to (p we get 



[a 2 (v,°'-»7 o )(-fl) 1/a ] J a=0. 



(23) 



Hence we see that the baryon current density nu a is con- 
served, 



K);„ = 0, 



(24) 



as it should be. 

Varying S with respect to a gives in turn (a prime 
signifies derivative with respect to a 2 ) 

a, a ' a ~a (<^, Q -V a )(^* -V«) - a V'(<j 2 ) = 0. (25) 

This last equations provides a bridge to the thermody- 
namic variables. As in the case of the superfluid, which 
emerges from complex scalar field dynamics only for long 
wavelengths, so here. We imagine that whereas the phase 
<p may be rapidly varying, the amplitude a is slowly vary- 
ing. More precisely, if L denotes the scale of a variation, 
we assume 



L- 2 «|(^-0(^ a -r? a )|. 



(26) 



But <T, a ' a /a is then Q(L 2 ), so it must be negligible 
compared to (ip, a —ri a )(ip, a — r) a ). We thus have 

(<p, a -ri a X<p, a -V a ) = -V'(o a ). (27) 

With the help of Eqs. l|2*2*|l and (fTTJl this can be written 

( 2 n 2 =a¥(a 2 ). (28) 

from which it is clear that we must require V' > 0; re- 
gions of spacetime where this is not true cannot be occu- 
pied by the fluid. 

C. Electric current and Lorentz force 



If we vary S with respect to A a we obtain 



f^ F «/3_^ u [« 6 /3]\ ( _ 5) l/2 
V47T 27T / 



= 0, 



or more conveniently 

(f q/3 -2u[ q 6 /31 ) =0. 



(29) 



(30) 



Comparing this with Maxwell's equations F a P -p = 
Att J a , where J a is the electric 4-current density, gives 
us 



j« = _L( u [«6ffl) 



(31) 



Carrying out the derivatives on suitably factored quanti- 
ties and taking cognizance of Eq. I|24|) we have 

4nJ a = (nu a ). g (—) + (—) nu a -(-) nvP . 



■3 



(32) 



Finally in view of Eq. Q we have for the Lorentz force 
density 



"-Fa/3 



n / ;7 
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(33) 



It is perhaps significant that the quantity in the square 
brackets is the Lie derivative of b a /n along n«°. 

With an eye on reconstructing the Euler equation from 
our approach, let us use this result to write the Lorentz 
force in an alternative form. We first define 



(34) 



We now substitute Eq. (|18|) for j] a here and carry out the 
derivatives: 

b 1 /6 7 \ 
47rC^« = (F art p-F/3 r> a) — \-F aj ( — J -Fpy I — J . 

"(35) 

By the Maxwell equations 

F a p- n + F ia - t fj + Fp TlCe = 



(36) 



we may replace F ai -p — Fp l;a in the previous equation by 
F a p- n . This done, let us interchange the a and (3 indices, 
and contract the resulting expression with nu^ . After 
completing a derivative we have 



4tt C T a a nuP 



(F a/3 mi?) — + F af3 (m/) ;7 — 



nu 



- -F 10 (-) 



.(37) 



By virtue of the MHD condition the first and fourth 
terms on the r.h.s. of this last equation drop out. By 
comparing with Eq. (|33|l we see that 



F a p J p = C F a p nu 
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(38) 



Armed with this identity we turn to the recovery of Eu- 
ler's MHD equation. 



D. The MHD Euler equation 

Let us introduce the notation 

v = (V') 1/2 , (39) 

where we explicitly mean the positive root. In view of 
Eqs. l|2*5)l and (|2*2*|l we now have 



n/u = C a , 



(40) 
(41) 



We now take the gradient of Eq. I|27|l and make use of 



vu a ((p. a p - rjc-p) = -w,p ■ 



(42) 
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We would here like to interchange the a and j3 deriva- 
tives of (p; however, we take note that the tensor 



W a p = tf-pa - cp ;a p 



(43) 



need not vanish identically, as we might suppose, because 
tp is a modular variable, and hence can be singular on 
various curves. We may thus rewrite Eq. as [see 
Eq. EH] 

V U a ((f.fja - W a g - 1]p. a + T a p) = -V V,{3 . (44) 

Dividing through by v, using l|41[) again, and observing 
that u a up ;a = Dufj/dr where r is the proper time of the 
observer with 4- velocity u a , we get 

—r J± = --U/3U a 1 - , (45) 

GtT V V V 

which has the form of a fluid's equation of motion. This 
has to be compared with the MHD Euler equation which 
has the form 



Pup 
dr 



P,0 a P,a . Fp a J a 

U U „ "I • (46) 



P + P 



p + p p + p 



First we shall work in a region outside singular curves 
of if where W a p = 0. The non-magnetic force terms on 
the r.h.s. of Euler's equation all contain the pressure 
gradient in a particular way. We observe that Eq. 145(1 
has terms with the same structure involving the gradient 
of v. This suggest that the two sets are identical. But 
this can be true only if 



dp dv 



p + p v 



(47) 



with the thermodynamic differentials taken at constant 
entropy. On the other hand, from thermodynamics we 
know that 



dp _ p + p 
dn n 



(48) 



where n, the baryon proper number density, is the same 
quantity we have been using above. Using the last two 
equations we have 

d(p + p) = p + P ^Ap_dv_ = p + p + p + p dv 
dn n dv dn n is dn 

Dividing through by (p + p) j dn reduces this to 

d{p + p)/(p + p)=dn/n + dv/v, (50) 
with integral 

v = P(p + p)/n = Pp. (51) 

(P is an obviously positive constant of integration). Thus 
if our guess is correct, v is proportional to the relativistic 
enthalpy per baryon of Eq. Q . 



Of course the T a $ term in Eq. (|4*5T> should be equiva- 
lent to the F al 3 term in Eq. l|4^|) . In light of result ijST]) 
this is consistent with the identity (|38|l provided 



p = c 1 . 



(52) 



We shall thus henceforth write p, — C,v. In view of the 
above we may conclude that in the region outside singular- 
curves, Eq. (|45|l from the field formalism coincides with 
the MHD Euler equation. 

Does this conclusion hold on the singular curves ? That 
would entail establishing that the constraint 



W a(3 u fJ = 



(53) 



holds on such curves. We shall indeed prove this last 
condition in Sec. IV Bl below. But for now let us note 
that formally Eq. H45J) is an integrability condition for 
Eq. (|41|) . First we notice that from the definitions of 
acceleration and Tp a , 



Du 



dr 



h Ufj U 

V V V 

_ „ [(^ U [a);/3] +??[q;/3]]w /3 



(54) 
(55) 



But according to Eq. l|4*T)> this is equal to —W a pu a jv. 
Hence Eq. I)45|l is satisfied identically. Actually this is not 
surprising; the method used here in deriving it involves 
"taking the curl" of the expression for u a , Eq. l|4*T)l . 



E. Gauge symmetry in the theory 

The vanishing of the divergence (|30(l means the quan- 
tity in brackets is the dual of a 4-curl. More precisely, 



pa _ 2u [a b l3] = */ a/3 , 



(56) 



where the antisymmetric tensor f a p and its dual *f a P are 
defined in terms of a field a a : 



fa(3 = 0-0, a — O-a,0, 

y a = i £ a ^ s f lS = e a ^ 5 a s ,^. 



Let us take the dual of Eq. (|56|l : 

*F<x0 — £ a 0jS u 1 b s = —faf} 



(57) 
(58) 

(59) 



As with any antisymmetric tensor, we can relate f a p to 
electric and magnetic-like vectors defined with respect to 
it' 3 , namely e a and b a given by 



ea = fa0U^, 

b a = -*f a pu fJ . 



(60) 
(61) 



Contracting Eqs. (|56|l with up shows that b a , as just 
defined, is identical with the Lagrange multiplier b a . And 
contracting Eq. i|5*§|) with u 13 shows that 



-a ■'-'Of 



(62) 
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We thus find that whereas F a p comprises only a magnetic 
part, f a p comprises both electric and magnetic parts, the 
first copying the magnetic part of F a p . 

The theory as stated has a U(l) x U (1) symmetry. The 
action and equations are invariant under 



Integrating this, and then integrating by parts leads to 



A a -> A a + X^\ a , 

&a ~> 0-ot ~f~ .(xi 



(63) 
(64) 



with A*- 1 ) and A^ 2 -* two independent functions. This sym- 
metry survives in the equations of motion after applica- 
tion of the constraints. 

We observe that the action S$ is invariant under the 
transformation 



tp- 



A, 



(65) 
(66) 



for arbitrary A. However, this gauge invariance is explic- 
itly broken by the 3c, part of the action. It is true that 
once one takes the conservation law (|24[) into account, 
any change of induced by transformation (|66|) can be 
converted into a surface term by means of Gauss' theo- 
rem. But the use of equations of motion derived from 
the action is inappropriate at a stage where symmetries 
of the action are being considered. Thus the U(l) x U(l) 
symmetry is all there is. 



IV. CHARACTER OF THE MHD FLUID 

There are two ways to bring out the character of the 
MHD fluid in the theory. One is to exploit Eq. (|TJ and its 
logical predecessors, all of which follow from differentia- 
tion of Eq. (2ZJ. This will be done in Sec. llVAl The sec- 
ond way is to formally construct the energy-momentum 
tensor as it would occur in Einstein's equations. This 
road will be travelled in Sec. IIV Bl 



A. Equation of state 

As with any scalar field representation of a fluid, we 
may establish the relation between the form of V(a 2 ) and 
the equation of state. In light of Eqs. (|39|l and Q52[l . we 
first write Eq. I|48|) in the form 



dp 



dn 



dn 



Then from Eq. (|4(J|I we obtain 
dn 



c^r = [v> 2 )] 1/2 + vvv)[vV)r 1/a 

dcr 2 



so that (|67|l becomes 



dp 
aW 2 



-a 2 V"(a 2 ) 



(67) 



(68) 



(69) 



p=-\V(a 2 )+a 2 V(a 2 )} 



(70) 



We now need an analogous expression for p. Substi- 
tuting Eq. (JTJ in Eq. (|47|l allows us to write 



dp 



dv 



da 2 ~ Cn da 2 



(71) 



Recalling that v = (V) 1 / 2 , carrying out the differentia- 
tion, and substituting here from Eq. I|40(l we have 

|£ = VVV). (72) 
Integration followed by integration by parts leads to 

p=-[a 2 V / (a 2 )~V(a 2 )}. (73) 

2 

Together with Eq. (j7U|l this equation provides the equa- 
tion of state of the fluid in parametric form. What are 
the general requirements on V(a 2 ) which yield a physi- 
cally acceptable equation of state ? One emerges from 
the need to enforce causality, specifically that the squared 
speed of purely acoustic perturbations, t> s 2 , be positive 
and less than unity. Dividing Eq. i|72|) by (|69|) gives 



dp 
dp 



V" 



V" + 2V'/a 2 



(74) 



Since we have already required V > 0, we must have 
that V" > as well. These two conditions automatically 
insure that p > 0. Requiring p > is perhaps too strong 
a condition: there is also electromagnetic pressure which 
may counter instability from negative fluid pressure. But 
it is safe to require that the strong energy condition be 
satisfied: p + 3p > 0. From expressions l|70() and l|73() 
this gives V < 2a 2 V, i.e., the logarithmic slope of Via 2 ) 
must be larger than 1/2. Thus the potential 



V(a 2 



A 



2 k 



(75) 



with nonnegative coefficients A 2 k and finite K, satisfies 
these conditions. 

Example 1: Assume V = (m/() 2 o 2 (m constant). 
This corresponds to p = 0, p = m and p = mn, and rep- 
resents a pressureless gas whose baryons have rest mass 
m. 

Example 2: Take V = /3a 4 . This gives p — -| p while 
p = 3 • 2- 5 / 3 C 4/3 /3 1/3 n 4/3 and fi = (y/Wa. This repre- 
sents thermal radiation with a sprinkling of baryons: the 
proportionality p oc n 4 / 3 is typical of adiabatic compres- 
sion of radiation. 



B. Energy- momentum tensor 

Einstein's equations are derived by varying the total 
action (including the Einstein-Hilbert part) with respect 
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to g a " . From this follows that the energy-momentum 
tensor T a p that sources these equations is determined by 

5S = - 1 -[T al3 (-g) 1 / 2 5g a P d 4 x. (76) 

where the variation is to be carried out before any of 
the constraints or equations of motion found earlier are 
enforced. 

Let us first carry out the variation of S c . Because 
A a ,b a and u a are regarded as fundamental variables, 
the only contribution of the first term comes from the 
g a P dependence of (— g) 1 ^ 2 - This will give rise to a g a p 
factor which will multiply Fjsb s u' y . However, once we 
take Eq. Q into account this term vanishes. The second 
term in S c depends on the metric in two ways; however, 
it is all multiplied by k which we know will vanish. Hence 
S c contributes nothing to the energy-momentum tensor. 
Strictly speaking the above is true only in the region oc- 
cupied by the fluid, that where n ^ 0. For without this 
condition we are unable to derive the MHD condition, 
or the condition n — 0. So we should say that S c con- 
tributes nothing to the energy-momentum tensor of the 
fluid. A similar remark will apply below. 

Variation of S nor with respect to g a @ produces two 
types of terms, both of which are multiplied by £, which 
we know will vanish. Hence S nor also does not contribute 
to the energy-momentum tensor of the fluid. 

Turning to we meet a similar situation. u a and rip 
are fundamental variables, so no metric enters in their 
scalar product. Variation of g a p in (— g) 1 ^ 2 contributes a 
factor g a P(—g) 1 l 2 which is multiplied by the u 7 ?y 7 . Since 
this last term ultimately vanishes by Eq. I|19|l . too 
makes no contribution to the fluid's energy-momentum 
tensor. 

It is plain that variation of Sm with respect to the 
metric will produce the Maxwell energy-momentum ten- 
sor. Therefore, the MHD fluid energy-momentum tensor 
comes exclusively from S^. Carrying out the variation 
gives us the sum T^O + T^f) , where 

T^p = a, a a, l3 ~-a, J a,' 1 g a p, (77) 

ee <r 2 (<p, a -v a )(.<P,(i-rifi) (78) 
- i [a 2 ^ n -77 7 )(^ 7 -V 1 ) + V(a 2 )} g a0 . 

Let us look at . The first term is 0{a 2 L~ 2 ) in the 
terminology of Sec. Ill CI According to inequality (|2*uTl . 

that term is negligible compared to the first term in T^) . 
By the same logic <r, 7 <r, 7 may be neglected compared 
to a 2 (ip n — ?7 7 )(<£, 7 — -q 1 ). We thus see that the fluid's 

energy-momentum tensor is dominated by Xj^ % . Let us 
replace in this last cp, a —r/ a by its expression (|41|l and 
simplify by means of Eqs. (|39[) - l|40l) . The result is 

Tap = °* ^V) "a U0 + - [a 2 V'(a 2 )-V(a 2 )]g aB . (79) 



This may be compared with the standard perfect fluid 
energy-momentum tensor 

T S = {p + PWu p +pg a p. (80) 

The p and p may now be identified; they coincide with 
those given by Eqs. 1)70(1 and l|7H|l. The pictures obtained 
from the field equations and from the energy momentum 
tensor arc thus consistent. 



V. THE 1+3 VIEWPOINT 

A. Preliminaries 

As we have seen, our proposed field theory for MHD 
leads naturally to the fluid 4-vclocity 

u a =(p a -r, a )/v, (81) 

and it gives rise to the projection tensor h a g = g a Q + 
u a U/3 projecting orthogonally to u a . Decomposing the 
covariant equations in parts parallel and normal to u a 
will render the effects of curvature more transparent and 
enhance our physical understanding. We adopt in the 
following the notation of Ref. Il3| for convenience. 

There are now two differentiation operators, namely 
the (proper) time derivative T> = u a W a (often denoted 
by an overdot: Vf = f) and the (totally projected) spa- 
tial derivative D a (e.g., D a Tg 7 = h a<r h\ V CT T/), 
respectively. The covariant derivative of any scalar func- 
tion / is thus decomposed as V a / = D a f — u a T>f. 
A congruence of observers with 4-velocity u a is covari- 
antly characterized in terms of the associated Raychaud- 
huri |l4j kinematical variables, which are derived from 
the covariant derivative of u a : 

V a «/j = -u a Uf3 + T> a up , (82) 
D Q u/3 = \Qh a p + <j a p + us a p . (83) 

Here the trace O = Y) a u a is the volume rate of expan- 
sion of the congruence; o a p = D( a up) — ^Qh a p is the 
trace- free symmetric rate of shear tensor {o a p = crr a a) , 
VaBU 13 — 0, a a a = 0), describing the rate of distor- 
tion of the congruence; and Lo a p = D[ a um is the skew- 
symmetric vorticity tensor (u) a p — u)\ a p], uiapu 13 = 0), 
describing the rotation of the congruence relative to a 
non-rotating (Fermi- Walker propagated) frame. Finally 
T>u a — ii a — u^V pu a is the relativistic acceleration vec- 
tor, which represents the influence of forces other than 
gravity on the observer (a free-falling observer has van- 
ishing acceleration in her rest-frame). 

A crucial feature of the newly introduced derivative 
operators is that they do not commute with each other 
in general. In particular, one finds for any scalar function 
/ the commutation relations 

D [a Dfl/ = (Jeff, (84) 
Vj-hJVppf) - —u a f + (DaU 13 ) Dpf . (85) 
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Eq. H84|) immediately tells us that the spatial derivative 
D is a covariant derivative if and only if the vorticity io a p 
vanishes, in which case the 4-velocity u a is hypersurface- 
orthogonal and h a p becomes the induced metric of the 
hypersurface. 

It is convenient to introduce the totally antisymmetric 
tensor (spatial volume element) e Q /3 7 = u s esa(3^ ■ With its 
help we define the spatial curl of a 3- vector V a (V a u a = 
0) to be cur\V a = e^^DpVy. Moreover, we define the 
vorticity vector ui a = ie"^ 7 w / 3 7 such that uo a p — e a p 7 uj' y . 
Analogously, B a = \t a P 1 Fp 1 and F a p = e a/ 3 7 i? 7 . 

We note that the above commutation relations have 
to be modified for the case of the modular field <p, since 
its curl W a p is not vanishing. Decomposing W a p into 
electric and magnetic components, 



W a/ 3 =u a Wp- up W a + e^V 1 , 



(86) 

the modified commutation relations l|84|) and l|85|) can 
then be written as 



D [a D /3]</? = £a/3 7 <P + \ V ~') 



(87) 



and 



T> a ip-hJV(T)pip) = -u a ip+(D a U )D ff ip+W a , (88) 

respectively. Note that in this section we do not en- 
force the constraint ((53(1 . W a = W a ^up = 0, in order to 
demonstrate which findings hold irrespective of it. 



B. Kinematical variables 

We now determine the kinematical variables associated 
with a congruence of observers with 4-velocity u a given 
by 1)81(1. First by contracting (|81J) with u a and h a/ 3, re- 
spectively, we get 



v = -if, 



Likewise, from Eq. (|4(J|I we obtain 
n = — C _1 <j2 ^P ■ 



(89) 
(90) 

(91) 



Hence the number density n is determined in terms of 
both a and (p, while the field r\ a is determined by the 
spatial dependence of ip alone. 

The expansion of the congruence is most easily ob- 
tained by writing out the divergence in Eq. I(24|) and re- 
alizing that D a u a — V a u a : 



e = -2. 



(92) 



The congruence u a has non-vanishing vorticity. This 
is most readily inferred by making use of the modified 
commutator relation (|87|l for p and also Eq. (|90|l . giving 



U a (3 = \ (^>[ a VP] - 7j £a/37^ 7 



(93) 



Further, the shear of the congruence is given by the ex- 
pression 

1 In 

<J a f3 = — rD( a [V^)^-))/3)] +--h a p, (94) 

ip y L ' 6 n 

where we used Eq. I(92|) to replace the expansion term. 

It remains to calculate the 4- acceleration u a . We use 
the Euler equation (|45|l for this task, which we write in 
the form 

i>u a = -D a is+((ny 1 F a p J l3 -W a , (95) 

where we have employed Eq. Writing F a p — 

CapjB 1 and replacing v via expression 1(89(1 . one gets fi- 
nally 
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B a p> + -^-e al3l JpB 1 -W a 



(96) 



It is instructive to calculate u a in another way, namely 
by employing the commutation relation 1(88|) with 19011 . 
Doing so one obtains 



- (D a ip - \ri {a) +n D a u 13 
p \ L 



W a 



(97) 



where the angle brackets denote projection with 
h a f3. Comparing the last two expressions for the 4- 
acceleration, we find an evolution equation for the field 
r) a , namely 

ri (a) + \<dr, a + (a a p + uo a p ) rf = - {Qn)" 1 e a ^J^B y . 

(98) 

We thus see that the electromagnetic interaction induces 
spatial inhomogeneity in the scalar field p via the Lorentz 
force J x B. The last equation will be useful in Sec. lVIII 
It may be worthwhile to investigate the relationship 
between the tensors T a p — 2 V[„ D^j ^ and W a p = 
2V[ Q ,V^]^, respectively. Expanding D^ in the defini- 
tion of T a (3 in terms of V/3 we readily obtain 



(99) 



This can be further rewritten by employing in the r.h.s. 
the decompositions (|5fj|l . the commutator relations 
(|57jl . (|8~8"|) . and the result (|9l5)l. Some straightforward 
algebra reveals that 



T a p = —u [a F /3]7 J 7 + 2 D [Q jjfl, (100) 
s n 

from which the electric and magnetic parts of J- a p are 
manifest. The magnetic part may be written in an al- 
ternative way by means of Eq. 1)87(1. Eq. ((100(1 is fully 
consistent with Eq. ((38(1 . 

We are now in the position to give a derivation of 
the constraint ((53(1 . that is W a = 0, by looking at 
overall energy-momentum conservation. We recall from 
Sec. IIVBI that the action for the complex scalar field ip 
yields a perfect fluid energy-momentum tensor T^J [cf. 
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Eq. I|78(l ] in the long-wavelength limit, with the number 
density n and the potential V related via Eq. (|2*S1) . The 
concomitant energy density and pressure of the fluid as 
given in (|7U)l and (|73^) can be recast as 



a 



Via 2 



(a^) 2 + V(a 2 ) ,(101) 



(a if) -V{a 2 ) ,(102) 



where Eqs. and (f^TTl) have been used. 

Now, the overall conservation of energy-momentum 
requires V pT^i — F a/3 Jp. The familiar energy and 
momentum conservation equations follow by contracting 
this with u a or by spatially projecting on the index a, 
respectively: 

p + Q(p + p) = 0, (103) 
(p+p)u a +D a p = e a ^JpB y . (104) 

We can check the consistency of these equations by in- 
serting the expressions H101(l and (|lt)2(l . Observing that 
V = 2aaV' (with an analogous relation for D a V) and 
remembering Eq. (|92|l . the energy equation 11031) is read- 
ily seen to hold. Now, by calculating the spatial gradient 
of the pressure, D"p, taking the commutator relation 1|85|) 
as well as Eq. <|9"S1) into account, and comparing the result 
with the momentum equation l|104|l . we end up with 



(nW a = 



(105) 



This means that, in regions occupied by the fluid, the 
constraint l|53|) must hold in order to implement the mod- 
ularity of the scalar field tp in a consistent manner. Thus 
we have shown that the derivation of the MHD Eulcr 
equation from Eq. J57J), as carried out in Scc. lIII Dl is also 
valid on singular curves. And, of course, the momentum 
conservation Eq. (|104|) is already the Euler equation, here 
obtained directly from energy-momentum conservation. 



C. The role of the Lagrange multiplier b a 

In order to get a better understanding of the "little 
magnetic field" b a , we look at Maxwell's equations from 
the 1+3 viewpoint. Due to our MHD condition of van- 
ishing electric fields, E a — 0, these become now 



4tt J a = F^.p = (e^B^.^ 
= *F af5 . l3 = -2 U a B® 



;0 



(106) 
(107) 



As usual, the 4-current is decomposed relative to u a in 
the manner 

J a = gu a + J {a) ; = h a p J 13 , (108) 

where g — ~J a u a is the charge density and is the 
Maxwell 3-current, respectively. Projecting the above 



Maxwell equations perpendicularly and along u a , one ar- 
rives at the following system of equations: 



4tt 

4-7T g 





curl r + e"" up J3 7 , 
-2uj a B a , 



(109) 
(110) 



3 6B Q -(a%- 0/^)^,(111) 
D„ B a . (112) 



4tt J q = F at3 . p = 



=*F' 



ap 



These equations are just the familiar special relativistic 
Maxwell equations with vanishing electric field E a but 
with the general relativistic correction terms encoded in 
the kinematical variables of the observers' 4- velocity u a 
(c.f. Ref. Ila. for example). On the other hand, in light 
of Eqs. (|56|) and H59|) Maxwell's equations may equally 
well be written in the form 

2(V Q 6 /31 ) ,(H3) 
( e "/37 5 7 _ f«py^ = _ 2 ( U I« e /3l) ^ .(114) 

It is immediately clear from Eq. I|113fl that the Lagrange 
multiplier b a plays a role analogous to an electric field, 
whilst Eq. i|114|) is identically to Eq. (|107fl since we al- 
ready know that e Q = B a . Decomposing the inhomoge- 
neous equation (|113f) readily yields 

AttJ^ = -6 <a> -|e& a + ( ( T Q /3 -w Q /3 )6 /3 ,(115) 
Ai:g = D Q b a , (116) 

which look like the Ampere-Gauss equations, and once 
again demonstrate the "electric" nature of the La- 
grange multiplier b a . Thus we might view (nr] a — 
(47r)~ 1 e Q/?7 bp B 1 as a Poynting vector. 

It should be noted that the Lagrange multiplier b a is 
anything but independent. Suppose suitable scalar fields 
a and <p were chosen, fixing u a and the kinematical vari- 
ables, and a solution for the magnetic field B a had been 
found. Comparison of Eqs. HlU9f) and (|115|) then reveals 
that b a can be determined via a differential equation 
wherein B a acts as a source. Our model is thus con- 
sistent as long as that differential equation possesses a 
non-trivial solution for b a . Such a solution is not unique. 
For according to Eq. (|lllfl . a multiple of B a may be 
added to b a without changing the latter's status as a so- 
lution of Eq. i|115|) . The addition to b a of a multiple of 
u a , which is suggested as a possibility by Eq. (|113p . is 
made untenable by Eq. (|12(l . 



VI. CIRCULATIONS AND HELICITIES 

A. Conserved circulation 

We return to Eq. (|81|l . Solving for ip <a , substituting 
rj a from Eq. I|18|) , and integrating the result over a closed 
loop initially lying on a spacelike surface and carried 
along by the flow, we obtain 



{pu a + CVa)dx a = 2tt(N, 



(117) 
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where N S Z counts the number of times the phase i/5 
winds around its natural interval [0, 2tt] as the point tra- 
verses the loop once. At a purely classical level this result 
gives conservation of the circulation To ■ This is because 
hydrodynamic evolution amounts to a continuous defor- 
mation of the loop and a continuous variation of the in- 
tegrand in To- By continuity ./V cannot change during 
such an evolution, and so To is conserved. 

The conservation of To was put in evidence by Beken- 
stein and A. Oron ^lj; they employed a purely hydro- 
dynamic variational principle for MHD (in contrast to 
the present field-theoretic one) to generalize a circula- 
tion conservation law discovered by E. Oron and given in 
Ref. (henceforth BEO) for flow with both stationary 
and axial symmetry. What is new in our result here is the 
quantization of the circulation, a direct result of the use 
of a phase as a dynamical variable. By contrast, Ref. ITTI 
used only thermodynamic variables and the velocity to 
describe the fluid. 

What is the meaning of the quantization ? How does it 
square with the presumed continuity in allowed values of 
circulation ? We notice from Eq. (|10f> that the dimensions 
of £ must be those of action. This because r] a carries 
the same dimension as 93, Q , namely those of reciprocal 
length, L^ 1 , while n by definition has dimensions of L~ 3 , 
and u a is dimensionless (c = 1). Now let us consider a 
situation with arbitrarily weak magnetic field frozen into 
a nonrelativistic fluid. According to Eq. the second 
condition means /1 equals the rest mass m per particle; 
and the spatial part of u a reduces to the usual 3-velocity 
v. The condition 1117|) is thus 

<pv-dt = N—. (118) 
J m 

Now our action (JJJ in the limit of vanishing magnetic 
terms, and with neglect of derivatives of a, is suitable for 
describing a superfluid condensate. However, for a super- 
fluid circulation is quantized by the Onsager-Feynman 
rule 0|, which is precisely l|118|) with £ = H. We must 
thus calibrate £ to a value h in all cases (although in def- 
erence to our classical approach we shall continue to use 
the notation £). The smallness of the quantum of action 
means that for macroscopic loops the Oron circulation 
To will take on an almost continuous range of values, as 
would be expected classically. 

The above approach is germane to a field theoretic 
approach; but how can one understand the conservation 
of To at the macroscopic level ? For this purpose we 
shall define the 4-vector 

w a = vu a + 7y Q , (119) 

this being nothing else than £ _1 times the vector consti- 
tuting the integrand of Eq. (|117ll . The fact that w a = ip^ a 
does not make the vector's curl W a p [defined by Eq. (@2J)] 
identically zero because, as mentioned earlier, (p is a mod- 
ular variable, and hence can be singular on various curves 
(these would be the cores of vortices of the spatial vector 



h a @ Wfj = rj a ). As shown by Eq. 

W af3 vP = 0. (120) 

An obvious identity satisfied by the curl W a /3 is [see 
Sec.HI 

WWry]=0 or W a % = 0. (121) 

We now consider the rate of change of Oron circulation, 
or more precisely, its change when the contour is Lie- 
dragged along u a . As a first step we write by means of 
Stokes' theorem 

r -c j>W aP dZ afi , (122) 

with d£ a ^ the 2-area element tensor on a surface span- 
ning the loop in question. Secondly, according to a cal- 
culation by Bekenstein and E. Oron (BEO) [12||. for any 
tensor W a p and any loop 

Cu T = cf [2(^ W l[(3 ) :a] - 3W [a0n] u>] d£ Q/5 . (123) 

And of course, the integrand here vanishes by Eqs. I|120|) 
and (|121fl . Thus the circulation To is conserved quite 
apart from the pertinent vector being the gradient of a 
phase. Both because of this and because of its limiting 
form when F a p — > 0, we may regard To as the general- 
ization of Kelvin's circulation to MHD. 

We recall that the magnetic (Alfven) circulation 

T A = <j> A a dx a = J F af3 d^ a(3 (124) 

is also conserved in perfect MHD. This may be shown 
exclusively from Eq. I|123|) with W a p — > F a p and the 
conditions 0} and (J2SJ). In Sec. IVII B"3l we shall exhibit 
a third, new, circulation which is conserved in situations 
with two spacetime symmetries. 

B. Conserved helicities 

According to Moffatt [lj| the conservation of fluid he- 
licity in perfect pure fluid dynamics reflects the conser- 
vation of Kelvin circulation around linked vortex lines. 
Likewise, the conservation of the fluid-magnetic Wolt- 
jer helicity m MHD reflects the conservation of the 
Kelvin circulation around a vortex and the Alfven circu- 
lations around a magnetic flux line linked with the vortex. 
And Woltjer magnetic helicity reflects the conservation 
of Alfven circulation around linked flux lines. Having 
displayed the new circulation To, we may ask whether 
it furnishes new helicity conservation laws to replace the 
lost Moffatt helicity in the MHD regime. 

For a perfect pure relativistic fluid the Moffatt helicity 
conserved current is [jji |2(j 

H a = i £ a^s ^ U[s y n] ^ = _^ a . (125) 
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This may be compared with the conserved currents for 
Woltjer's magnetic and fluid-magnetic helicities in per- 
fect MHD flow [Ullj: 

H u = A p = -B l3 Af3U a -cj)B a , (126) 

Hf m = *F a/3 imp = -^B a , (127) 

where <fi = —A a u a and A a = hj 3 ''Ap. In MHD the cur- 
rent Hf is no longer conserved. We propose to replace it 
by [notation as in Eqs. 1119(1 and 

H? = *W a0 wp. (128) 

The proof that H" is conserved is simple. In view of 
Eq. (|121|) we have 

H?, a = \ *W ai3 W a(3 , (129) 

where the antisymmetric part of wp ;a is the only part 
that survives contraction with the antisymmetric *W a0 . 
Now W a = by Eq. {53), so from Eq. JSE) we see that 

W a p = ec^su^V 5 , (130) 

*W a0 = -UaVp+UpVa. (131) 

It is obvious that the full contraction of indices between 
these two tensor vanishes identically; hence by Eq. I|129|) 

Hf- a = 0. 

Two additional candidates for conserved helicity cur- 
rents suggest themselves. One is *F a @ wp] however, in 
view of Eqs. 10, and (J^, this is no different 

from Moffat's fluid-magnetic helicity i? f " r The second 
is *W a/3 Ap. Writing *W a p in terms of wp and *F a p in 
terms of Ap gives 

*W af3 Ap - *F af3 wp = \ (e a ^ 5 w s Ap) . (132) 

Since the r.h.s. is divergenceless (by antisymmetry of 
£ a0fS Ws Ap), the proposed helicity current is indeed con- 
served as a result of the conservation of H^. 

It may also be seen that the conserved helicity (the 
zeroth component of *W al3 Ap integrated over space) dif- 
fers from the Moffat's fluid-magnetic helicity (the time 
component of Hf m integrated over space) by 

i J ( e <*fhS Wg A p )^dE a = \ I e^wg ApdY, ai 

(133) 

where dY, a is a 3-volume element on the corresponding 
spacelike surface while c?S Q7 is a 2-area element on the 
latter's boundary at infinity. This looks gauge depen- 
dent. However, assume we are able to isolate a "physical" 
part of Ap. We shall now show that the corresponding 
integral vanishes if the fluid extends to infinity. 

It is convenient to think of the boundary at infinity 
as spherical. Thus cE Q7 comprises only temporal and 
radial components. By the antisymmetry of e a01& only 
the angular components of wg and Ap will contribute to 
the second integral in Eq. I|133|) . 



Now asymptotically £ a p~/S ~ r 2 ; accordingly the rele- 
vant components of dS Q7 vary as r 2 while e Q / 3 T 5 ~ r~ 2 . 
Thus the integral has the asymptotic behavior of wg Ap 
with P and S angular coordinates. Now the phase tp 
should depend on angular coordinates and not fall off 
with radius, so wg does not decay with r. However, Ap 
must fall off with increasing r. Were this untrue, F a p 
with both indices angular variables would survive at large 
r which would mean that the physical radial magnetic 
field falls off as r~ 2 or slower. But this would imply a 
nonzero magnetic monopolc field which we may discard 
as unphysical. Thus the integral in Eq. (|133fl vanishes. 
It follows that the current *W a @ Ap does not furnish a 
new conserved helicity. 

In conclusion, perfect MHD has three different con- 
served helicities corresponding to the currents defined in 
Eqs. (|126|) - (|128|l . The fluid helicity conservation law con- 
structed on the basis of Oron circulation was first claimed 
in Ref. 20i, but only for stationary axisymmetric flow. We 
have now lifted the symmetry restriction. By contrast the 
claim in Ref. |53 that there is a pair of conserved fluid- 
magnetic helicities based on Oron's circulation has been 
here nullified by our demonstration that these are iden- 
tical with the conserved Woltjer fluid-magnetic helicity. 



VII. BERNOULLI THEOREMS 



The Bernoulli theorem is well known from nonrelativis- 
tic perfect fluid mechanics; it actually appears in two 
types [IJ. The type-1 theorem obtains in any station- 
ary adiabatic flow: the sum of the specific enthalpy, spe- 
cific bulk kinetic energy and gravitational potential is in- 
variant along streamlines. The type-2 Bernoulli theorem 
states that for isentropic stationary potential flow, the 
above mentioned sum is a constant all over the flow. Rel- 
ativistic formulations of both types are known 0, [2^] . 
Less well appreciated is the existence of an additional 
Bernoulli theorem for each space symmetry shared by 
the flow and the gravitational potential (or geometry in 
the relativistic form). 

Relativistic MHD also exhibits Bernoulli theorems. 
For example, in any flow with both time and azimuthal 
Killing vectors, BEO exhibited a pair of type-1 Bernoulli 
theorems 01 • It i s quite clear from their work that 
similar theorems arise if the mentioned symmetries are 
replaced by some others, but that two symmetries are 
needed together. In the sequel we show, by exploiting 
the field description of this paper, that one can actually 
obtain one Bernoulli theorem of type-1 and one of type-2 
for each spacetime symmetry of the geometry and flow. 
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A. Case with one spacetime symmetry 

1. Type-1 Bernoulli theorems 

Assume that the metric and the fluid 4- velocity possess 
a symmetry described by a single Killing vector . This 
£ a satisfies Killing's equation 

V Q ^ + V /3 r=0. (134) 

Let us decompose the Killing vector as follows: 

C = ku a + k a (k = -Cu a ; k a u a = 0). (135) 

Inserting Eq. (|135fl into Eq. (|134fl and contracting in turn 
with u a up and with u a produces the independent 
equations 

k + k a ii a = 0, (136) 
kii a + k {a) ^B a k-k> 3 'D a u = 0, (137) 

which we shall employ presently. 

Further we record the Lie derivative of a scalar function 
/ and of the 4- velocity field u a with respect to £ Q : 

2tf = f = k.f + k a D a f, (138) 

j2£u Q EE ^V U° - U Vf3 C Q 

= - (k + k 13 iip) u a + k p Dp u a - k {a) . (139) 

Since the flow partakes in the symmetry £ a of the space- 
time, u a = 0, which upon using the first Killing equa- 
tion I136[l in (|139[> . implies 

£>> = k p Dpu a . (140) 

Eliminating k^ between this equation and Eq. I|137|) 
provides the convenient expression for the fluid's acceler- 
ation 

ii a =D Q \nk + 2k~ 1 uj afs kp, (141) 

which tells us that in the absence of vorticity the accel- 
eration is the gradient of the "potential" — ln(£ Q w a ). 

Let us compare this result with Euler's equation l|95|) 
with W a = and Eq. 

u a = — D Q In n + (n^)- 1 e Q/37 JpB- t . (142) 

We infer that 

(pl^e* JpB 1 = B a h-i(kfi) + 2k- 1 uj a ' 3 kf 3 . (143) 

Contracting this with kp gives 

n- 1 ^ 1 k a J (i B 1 = (i k a D a In (k fx) = -V (k fi) , (144) 

where Eq. (|138fl has been employed to achieve the second 
equality. The magnetic term here can be gotten from 
Eq. I|98|) contracted with k a . The expression involving 



kinematical variables appearing therein can be deduced 
by contracting Eq. I|83l) with k a rP and simplifying with 
use of Eq. lfTiU|) . The result is 

n~ X e a ^ k a Jp S 7 = -({v^ k a + V a k {a) ) 

= -V((„ a k a ). (145) 

with the last equality following because both k a and r\ a 
are orthogonal to u a so that f]^ can be replaced by T>r] a , 
etc. 

Combining Eqs. (|144fl and l|145|) we arrive at the 
sought result: 

V(C,rfk a - kfi) = 0, (146) 

which is the relativistic MHD Bernoulli theorem. In other 
words, we have shown that 

^{u a C) + {^n)" 1 e a ^ i a bpB 1 = A, (147) 

where K is constant along streamlines. It seems not to 
be possible to prove that A" is a global constant in the ab- 
sence of additional assumptions, so what we have found is 
a type-1 Benoulli theorem for each spacetime symmetry 

2. Type-2 Bernoulli theorems 

Now choose coordinates in such a way that one of them, 
x^ , increases along the integral lines of £ Q . As an analog 
of what is referred to as potential flow, we choose the 
following ansatz for the phase in our formalism: 

ip. = Ax« + H{x R ) . (148) 

Here H is some function of the coordinates other than x^ , 
which we denote collectively by x , and A is a constant. 
We now show that this choice of ip is consistent with the 
postulated symmetry. 

It is clear that ip^ a does not depend on x^. Thus be- 
cause u a shares in the symmetry, it follows from Eqs. I|89|) 
and (|51|l that v and fi are x^ independent. From Eq. I|91|) 
we find, likewise, that n/a 2 is x^ independent. We then 
gather from Eqs. I|101|) and l|102|) summed together that 
p + p shares in this symmetry and then Eq. shows 
that n and thus a 2 do likewise. Returning to Eqs. I|1U1|) 
and l|102|) we verify the symmetry of p and p separately. 
All the thermodynamic quantities thus share in the sym- 
metry. Finally from Eq. (|90|l follows that r/ a is x^ in- 
dependent. Eq. (|81|l then confirms that u a comes out 
x^ invariant as originally assumed. The above argument 
does not prove that the ansatz l|148|) is unique, but only 
that it does not conflict with expectations for the as- 
sumed symmetry. 

We now argue that the ansatz 1|148[) indeed represents 
potential flow. Evidently u a should be a single valued 
vector; this requires single- valuedness of H(x R ). If in 
addition the coordinate x^ is not compact, e.g. is a 
time translation Killing vector, then ip is single valued 
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and the flow it represents does not contain circulation; 
it is potential flow. Much the same can be said when 
x^ is compact, as when £ a is an axial symmetry. All 
loops around the axis of symmetry have then the same 
circulation, and it is still possible, with care, to use a 
(multivalued) potential for this velocity field just as one 
may use a multivalued scalar potential to describe the 
magnetic field around an electric current filament. 

Let us now substitute the ansatz (|148|l into the con- 
traction of Eq. 1(81(1 with £ Q : 



nu a e + tv a e = k. 



(149) 



Because K here is a global constant, this is the promised 
type-2 Bernoulli theorem. To recast it into an alterna- 
tive form more reminiscent of nonrelativistic versions of 
Bernoulli's theorem we substitute — u a from the 
above equation into u a u a = —1 to obtain 



g urus 



+ tl - 1 (K-(r ]a ?*)g 



(150) 



We observe from Eq. 118(1 that £n a £ a is ( independent as 
well as quadratic in the electromagnetic field [according 
to Eq. I|115|) , b a is linearly related to the magnetic field] . 

Let us pass to nonrelativistic MHD flow in a static, 
nearly flat spacetime. Focusing on the time transla- 
tion symmetry, £ Q = 8f, we may write the squared 
space velocity as v 2 = g RS URUs, g^ R = g tR ~ 0, 



r 



— (1 — 2$n) ($jv is the Newtonian poten- 
tial), and p = run + e [m is the rest mass per baryon and 
e the internal energy density). We see from Eq. 1149(1 
that in slow motion with negligible pressure and negli- 
gible magnetic effects, K ss —to so that we may define 
£ by K = — m(l + £), and assume \£\ <C 1. Similarly 
we should think of the nonrelativistic specific enthalpy 
h = (e +p)/mn and C 7 la^ a / m a s small compared to 
unity (c = 1). It then follows from Eq. 1(150(1 to first 
order that 



2 1 



+ h + $ w - (c/mto Q r 



(151) 



This is the generalization of the classic nonrelativistic 
Bernoulli theorem to MHD. We see that £ is the nonrel- 
ativistic energy per unit mass. 

In order that the above results be really useful we must 
calculate r] a £ a 1 or equivalently F^pb^ . But we have not 
succeeded in solving Eqs. ((115(1 and 1(109(1 for b a . This 
remains a problem for the future. However, in a situation 
with two symmetries, £° and we can reach useful 
conclusions by comparison with BEO's type-1 theorem. 
We now turn to this. 



B. Case with two spacetime symmetries 

1. Type-1 Bernoulli theorems 

We shall here recast BEO's results in a transparent 
form. When the geometry and flow display two indepen- 



dent Killing vectors, £f and £j j BEO find that 

F a p$ = AF af i$, (152) 

= F afi ^e 2 , (153) 
B a = -Cn{uf3X P u a +x a ) , (154) 
B 2 = C 2 n 2 [(u 0X P ) 2 +X a X a ], (155) 
where A and C are constants along streamlines, and 

x a = - A& . (156) 

(The dimensions of A in the definition of \ a depend on 
the type of symmetries involved.) 

BEO also find the Bernoulli theorems (i = 1, 2) 



B 2 u a + Cn upx 13 Ba 



Aim 



(157) 



where K\ and Ki are again constants on each streamline. 
BEO show that 



[iu aX a = K 1 -AK 2 . 



(158) 



In all this BEO assume that the field component 
F a p £2 vanishes asymptotically; Maxwell's equations 
then force it to vanish identically. 

Let us now compute Bp (J? in Eq. 1(1571) by means of 
Eq. ^ 154(1 . and then eliminate C 2 n 2 in favor of B 2 by 
means of Eq. I(155|) . After simplification and a cancella- 
tion we obtain the type-1 Bernoulli theorems (i = 1, 2) 



B 2 
Ann 



[UaX 



a\2 



XaX 



Ki. (159) 



2. Type-2 Bernoulli theorems 

Let us now imagine that the flow is described by the 
ansatz 



K lX ^ +K lX ^ +H(x R ), 



(160) 



where K\ and K 2 are constants, a;^ 1 and x^ 2 are the 
coordinates along integral lines of £f and j respectively, 
and the x R are the two other coordinates. By the same 
method leading to Eq. ((149(1 we now get the two type-2 
Bernoulli theorems [i = 1,2) 



(161) 



Obviously K\ and K 2 are also constant along stream- 
lines, and so these two laws may well correspond to the 
two theorems, Eqs. (1159(1 . which certainly remain valid 
under the restricted flow described by Eq. ((160(1 . (Note, 
however, that we cannot regard A or C as constant all 
over the flow.) 

To back up the above identification we immediately 
verify by use of the definition ((18(1 that, regardless of the 
form of b a , 



V a (i? = Ar, a & 



(162) 
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This is because 



(163) 



A similar expression applies with , and so Eqs. (|152|) - 
l|153jl give Eq. 1(162(1 for any b a . If the identification of 
Eq. i|161|) with BEO's result is justified, we must find the 
ratio of the magnetic term in theorem (|159fl with i = 1 
to that with i = 2 to be exactly A, just as predicted by 
Eq. 1(161(1 in light of Eq. 1162(1 . This is easily verified if 
cognizance is taken of Eq. (|156fl . Finally, if we substract 
A times Eq. (|161fl as applicable with from the same 
equation for £f , Eq. (|158|l of BEO emerges. Thus, the 
identification with BEO is consistent, and the conserved 
expressions in Eq. I|161|) and 1159|) must be the same. 
Accordingly, we find in our formalism that 



B 2 

47m 



XdX P ~ U a X a Xp€j 
(UaX a ) 2 +XaX a 



(164) 



which exhibits the expected quadratic character of r\ a in 
the magnetic field. Putting this in Eq. I|161|) we get two 
type-2 Bernoulli theorems in explicit form. 

Example 1: Consider a situation with both a time 
Killing vector, £f = <5", and a translational spatial 
Killing vector in the x direction, = We shall 
work nonrelativistically from Eq. (|151fl interpreted as a 
type-2 theorem, and assume that the magnetic energy 
per baryon, i? 2 /87m, is not large compared to the ki- 
netic energy per baryon v 2 /2. Then it is unnecessary to 
take into account terms of 0(v 2 ) in the square brack- 
ets in Eq. 1164(1 which would generate terms compara- 
ble to 0(v 4 ), the like of which have already been ne- 
glected in Eq. (|151|l . By the same token we drop in 
Eq. (1 16411 corrections to the Minkowski metric which are 
of 0($n) — 0(v 2 ). The resulting Bernoulli theorem to 
0{v 3 ) is 



B 2 



N 



A 



Aimm V %Vx + A 



£. 



(165) 



We recall that in the last few equations A is constant 
only along streamlines. 

Example 2: Still starting from Eq. I|151fl we shall re- 
place the above £f by an azimuthal Killing vector (cylin- 
drical coordinates {r,z, ip}): = The streamline 
constant will here be denoted A. In the nonrelativis- 



tic limit we shall put u t k, —\ and 



Hi. 



where Q denotes the azimuthal angular velocity dip/dt. 
Again we neglect terms of 0(fl 2 ) in the square brackets 
in Eq. (|164fl . so obtaining 



b 2 



b 2 ( n + A 



N 



Annm \ 20 + A 



= £. 



(166) 



Example 3: We return to the situation with time and 
x symmetries. In Eq. <|161fl with £f = £f — we shall 
use the expression 1(164(1 . Nonrelativistically we approxi- 
mate it Q £" ~ — 1 and u Q £f ~ v x where v x is the compo- 
nent of the ordinary velocity in the symmetry direction. 



We take the metric everywhere as Minkowski's and ne- 
glect h in comparison with unity. The result is 



B 2 



= V, 



(167) 



47mm \A(2v x + A) 
where V = K 2 /m is the linear momentum per unit mass. 

3. New conserved circulation 

Eq. ((164(1 gives the projections of C,r) a onto the two 
Killing vectors. Can we reconstruct (^r] a fully from this ? 
For this purpose we could try finding the projections of 
rj a onto two independent vectors, both orthogonal to the 
plane spanned by £f and ■ We use the following two: 



^ 7 S^ U 



"f s u a - 2u /5 S / g 7 S Q7 



7 a/3 



€2 cf 



(168) 
(169) 
(170) 



It is immediately verified that Qf ^ Q = Q 2 £,i a — for 
i = 1,2. Now since Qf is a linear combination of u a ,£" 
and £2 i ^ follows immediately from the full antisymmetry 
of the Levi-Civitta tensor that Q" and Qf are orthogonal 
and hence independent as required. Therefore, it is no 
curtailment of generality to write 



Xf3X 







u p xp X° 



aiQI+azQZ, 



Ann [ (m 7 X 7 ) 2 + X1X 1 

(171) 

since one recovers from this ansatz both projec- 
tions 1(164(1 . and it has enough freedom left to represent 
any vector in the space complementary to the symmetry 
directions. Both scalars ct\ and a 2 must depend quadrat- 
ically on B a because rj a is quadratic in it. In fact, by 
demanding n a u a — and noticing that Q 2 u a — we 
find that 



01 



B 2 /{Aim) 



H a /3 £1 



a/3 



■2(u a E^) 2 



(172) 



We shall now argue that a 2 must vanish on grounds 
of parity. By Eq. 1(90(1 the space part of t] a must be 
a true 3-vector (changing sign under a space inversion) 
since ip should be a true scalar — invariant under inver- 
sion. Both u a, s space part and those of spatial Killing 
vectors must also be true 3- vectors. And by Eq. ((152(1 . 
A must be a true scalar. Hence any scalar product like 
w Q £f j u a X a , £,ia or \a X a must be a true scalar. Thus 
the square brackets in Eq. ((17111 enclose a 4- vector whose 
space part is a true 3-vctor. Likewise, the space part of 
Qi is a true 3-vector, while a\ is a true scalar. It follows 
that a 2 Q 2 must be a true 3-vector. However, it is clear 
from the second form in Eq. I(169|) that the space part 
of Q 2 is actually a pseudo 3-vector, so that a 2 should 
be a pseudoscalar. Yet no such can be built out of B a 
(whose space part is a pseudo 3-vector) and the £f which 
is also quadratic in the B field. The closest we come is 
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B a it ■ {^PisB^Ci^i) with i = 1 or i = 2. However, the 
factor in parenthesis here is just Fk/sCfC^i which van- 
ishes by Eq. (|153fl . We thus conclude that the Q2 term 
in (yf must be absent, and that we have reconstructed 
the full rf field. 

The above does not inform us further on Bernoulli the- 
orems because Qf is orthogonal to both Killing vectors. 
But somewhat surprisingly it gives a further circulation 
conservation law. According to Eq. (|117|) . the circulation 
of iJLU a + Qria is conserved, and according to BEO the cir- 
culation of ixu a plus just the first term on the r.h.s. of 
Eq. I|171(l is separately conserved. We conclude that 



ai Qi a dx a 



(173) 



is a new conserved circulation for MHD flows with two 
spacetime symmetries. 

Example: Consider again the situation with both a 
time Killing vector = 5f and a spatial one in the x 
direction = 6" . In flat spacetime but working fully 
relativistically we have S 7 ^ S 7<5 = —2 and M^S ( g 7 S Q7 = 

«t£l - Ux&- Thus Ql = - 2u " - 2u £l + 2Ux& SO 

that Q\ a dx a = —2(uydy + u z dz) where y and z are the 
usual Cartesian coordinates orthogonal to the symmetry 
direction. From the normalization u a u a = — 1 it follows 



that {u x - u t £?) 2 = u 2 - u\ = 1 + v 2 -\ 
these pieces it follows that the circulation 



£. With all 



is charged, it is represented by a complex scalar field 
minimally coupled to the electromagnetic gauge poten- 
tial. In this paper we have provided, in the framework 
of general relativity, a complex scalar field representation 
for the flow of highly conducting but neutral magnetized 
perfect fluid (perfect magnetohydrodynamics). The cou- 
pling to electromagnetism is via a vector field distinct 
from the electromagnetic gauge vector, but the theory 
has a full U(l) x U(l) gauge symmetry. The scalar field's 
self-interaction determines the fluid's equations of state. 

The principal advantage of this theory is that it leads 
directly to Oron's conserved MHD circulation which 
plays the role analogous to that of Kelvin's circulation 
in pure fluid flow. Here we obtain the Oron circulation 
without the need to assume two spacetime symmetries. 
Additionally, with two symmetries present we find an en- 
tirely new circulation involving the magnetic energy per 
particle and the velocity field. We have discussed the 
structure of the conserved helicity current which is asso- 
ciated with Oron's circulation. Finally we have put in 
evidence the existence of a pair of Bernoulli-like theo- 
rems which are associated with each type of spacetime 
symmetry exhibited by the MHD flow. In the case of 
two simultaneous symmetries we are able, by comparing 
with old results, to obtain explicit forms for the Bernoulli 
conserved quantities. 



B 2 u.ydy + u z dz 
Attn u 2 + u 2 



(174) 



around an arbitrary loop moving with the flow is con- 
served. 
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